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Abstract 

Canonical formulation of higher order theory of gravity requires to fix (in addition to the metric), the 
scalar curvature, which is acceleration in disguise, at the boundary. On the contrary, for the same purpose, 
Ostrogradskis or Diracs technique of constrained analysis, and Horowitz formalism, tacitly assume velocity (in 
addition to the co-ordinate) to be fixed at the end points. In the process when applied to gravity, Gibbons- 
Hawking-York term disappears. To remove such contradiction and to set different higher order theories on 
the same footing, we propose to fix acceleration at the endpoints/ boundary. However, such proposition is 
not compatible to Ostrogradskis or Diracs technique. Here, we have modified Horowitzs technique of using an 
auxiliary variable, to establish a one-to-one correspondence between different higher order theories. Although, 
the resulting Hamiltonian is related to the others under canonical transformation, we have proved that this 
is not true in general. We have also demonstrated how higher order terms can regulate the issue of branched 
Hamiltonian. 


1 Introduction 

We often encounter higher order theories, for example, in radiation reaction [1], noncommutative quantum field 
theories [2] , anyons (quasi particles appearing in two-dimensional systems, which are neither fermions nor bosons) 
[3] and particularly in the context of gravitation, leading to field equations with fourth order or even more deriva¬ 
tives. A quantum theory of gravity in any of its form, (eg., arising in the weak energy limit of a yet complete 
theory, may be the string theory, M-theory or supergravity) contains higher order curvature invariant terms 
leading to higher order theory of gravity. Higher order theories also appear in classical mechanical problems such 
as Pais-Uhlenbeck oscillator [4]. A host of canonical formalisms appear in the literature to handle higher order 
theory. Some of these lead to the same phase-space Hamiltonian, some different. Although, all these Hamiltonian 
lead to correct classical field equations, some of these are very different and not related to the others through 
canonical transformations (see appendix). Further, there happens to be a mismatch in different canonical for¬ 
malisms as far as boundary terms are concern. Any higher order theory is associated with boundary (end point) 
terms which do not vanish fixing only the co-ordinate at the end points. For fourth order theory, it is required to 
fix either the velocity or the acceleration at the end-points (boundary) in addition. However, it does not mismatch 
Cauchy data with the endpoint (boundary) data. This is because, canonical formulation of higher order theory 
requires additional degree of freedom. Ostrogradki’s technique [S], Dirac’s constrained analysis [B] and Horowitz’s 
formalism (originally developed for the theory of gravity) [7] fix velocity at the end points. 

In case of gravity, no nontrivial lagrangian Cg can be constructed from the metric g^^i, and its first derivatives 
alone. Even for Einstein-Hilbert action, Lagrangian depends on second derivatives of the metric, in the form 
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Cg(g,dg,d^g) ■ Under metric variation, setting = 0, Einstein-Hilbert action yields a surface term, so that 

the complete action reads 
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and the supplementary surface term is known as Gibbons-Hawking-York (GHY) boundary term [S]. In the above, 
£m, K and h are the matter Lagrangian density, the trace of the extrinsic curvature tensor and the determinant of 
the three-space metric respectively. Note that, one is not allowed to fix the extrinsic curvature tensor (Jfy) at the 
boundary to get rid of the GHY term. This is because, the correspondence between Gauchy initial value problem 
of field equations and the boundary value problem arising out of variational principle requires that the number 
of initial data should coincide with the number of boundary data. If Kij is fixed at the boundary, boundary 
data exceeds the number of configuration space variables, which mismatch Gauchy data. Ganonical formulation 
of general theory of relativity (GTR) has been performed successfully, taking into account the GHY term, and is 
known by the name of ADM formalism [^. 


Ganonical formulation of higher order theory of gravity requires additional degrees of freedom viz. (/ly, Kij ) 
nni, as already mentioned. Metric variation of higher order theory of gravity in its simplest form, viz., f{R) oc 
aR + /3i?", leads to a boundary term which vanishes provided Kij - the velocity in disguise, is kept fixed at the 
boundary. Mannheim m and Mannheim-Davidson m found no trouble regarding quantization holding position 
and velocity fixed at the end points, and apparently there is no problem as such. However, there are at least four 
reasons why, instead of Kij one should fix the scalar curvature R at the boundary, and supplement the action 
with a boundary term, S = 2/3 f Kf'{R)y/hd^x , in addition to the GHY term. Firstly, the above boundary term 
reproduces the expected ADM energy upon passing to the Hamiltonian formalism [13) . Secondly, fixing Kij at 
the boundary, one looses the most cherished GHY term, which is identified with the entropy of a black hole. 
There is practically no physical interpretation available, why the concept of entropy of a black hole should get 
lost in strong gravity, and reappear when gravity is weak, particularly, when the correct expression of entropy of a 
Schwarzschild black hole has been found in the semiclassical limit m- Next, f{R) gravity has been identified with 
scalar-tensor equivalent form in Jordan’s frame of reference, under redefinition of f'{R) = $ and R = , where 

prime denotes derivative with respect to the Ricci scalar R, or in Einstein’s frame of reference, under conformal 
transformation g^^i, = f'{R)giLv = , where the conformal factor to is related to an effective scalar field (j) by 

the relation w = Under variation, field equations are obtained keeping $ fixed in Jordan frame and ^ fixed 

in Einstein’s frame, at the end points. This is equivalent to fix R at the end points. This means scalar-tensor 
equivalence of higher order theory of gravity may be established, only under the condition, 6R\d\> = 0, which 
requires the supplementary boundary term, noted above. Finally, since Kij is the basic variable and also fixed at 
the boundary, it attains the same status as . So, there is no reason as to why the action should not be varied 
with respect to Kij too. Nevertheless, classical solution in the process, is restricted to just one possibility or even 
worse, to zero possibility M- Therefore, we suggest that instead of fixing Kij , the Ricci scalar R should be kept 
fixed at the boundary i.e. SR\dv = 0, and the action should be supplemented by an additional boundary term 
in the form, E = 2/3 f Kf'{R)'/hd^x . Now, since {hij,Kij) should be treated as basic variables for higher order 
theory of gravity, so in the Robertson-Walker metric (say for example) 


ds^ 


—dt -fa (t) 
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+ r^{d9^ + sin^ed(l)'^) 
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hij = 


R = 6 



z (say), Kij oc 2ad = z, 



(3) 


essentially the velocity (i) is treated as an additional basic variable, while (JRjav = 0) implies acceleration (z) 
is kept fixed at the boundary in addition to hij = z. Therefore, we propose that in order to treat all higher order 
theories on the same footing, acceleration should be kept fixed at the boundary instead of the velocity. In section 
2, we shall show in the context of Pais-Uhlenbeck action [4], how Ostrogradsi’s [5] and Dirac’s technique [6] fail 
to associate such proposition. Further, we show that under above proposition, Horowitz’s technique [7], is not 
able to take into account the boundary terms in general. In this context, we have modified Horowitz’s technique 
to construct a viable Hamiltonian formulation of a modified Pais-Uhlenbeck oscillator action and followed the 









procedure for higher order theory of gravity. 


Canonical formulation of higher degree theory is yet another problem, which has not been resolved uniquely 
as yet. If the Lagrangian contains velocities beyond quadratic, momentum does not appear linearly in velocities, 
resulting in multivalued Hamiltonian with cusps, usually called the branched Hamiltonian. This makes classical 
theory unpredictable, since at any time one can jump from one branch of the Hamiltonian to the other. It also does 
not allow standard canonical formulation of the theory. Further, as energy has to be an observable, the multivalued 
Hamiltonian does not allow momentum to ensure a complete set of commuting observable. Finally, it also suffers 
from the disease of having non-unitary time evolution of the quantum state. Related issues arise in cosmological 
models in extensions of Einstein gravity involving topological invariants and in theories of higher-curvature gravity 
m For example, higher order theories of gravity are often plagued with non-unitary time evolution in the weak 
energy limit when expanded perturbatively about the flat Minkowski background, although the general theory 
might be free from such disease. To get rid of this unpleasant situation, even if one constructs a theory with a 
particular combination (Gauss-Bonnet) of the curvature invariant terms, the resulting Lanczos-Lovelock theory of 
gravity m is plagued with the problem of branched Hamiltonian. Such unpleasant issue arising out of branched 
Hamiltonian was addressed long ago m- Starting from a toy model, 


^2 




dt, 


(4) 


Hennaux, Teitelboim and Zanelli m had shown that in the path integral formalism one can associate a perfectly 
smooth quantum theory which possesses a clear operator interpretation and a smooth, deterministic, classical limit. 
Nevertheless, it puts up question regarding the canonical quantization scheme. Further, such formalism can’t be 
extended in a more complicated theory beyond the toy model (|3]). The same issue has also been addressed by several 
other authors in the recent years [n muni HI]. However, in order to solve the problem, they had to tinker with 
some fundamental aspect viz., loosing the Heaviside function to obtain manifestly hermitian convolution together 
with the usual Heisenberg commutation relations |18j , sacrificing the Darboux coordinate to parametrize the phase 
space m, and the standard Legendre transformation [20) . Additionally, the Hamiltonians obtained for the same 
above toy model following the above two prescriptions [50] and m differ considerably, and not related through 
canonical transformation [22] . Therefore none of these techniques are fully developed or rigorous. However, one 
can bypass the problem by incorporating higher order term. For example, if action Q is associated with a higher 
order term in the form 


A, = 


cxf + 


dt, 


(5) 


the problem associated with branching is resolved uniquely, in the process of canonical formulation of higher order 
theory. This has been demonstrated recently in the context of Lanczos-Lovelock gravity [22] . 


In a nutshell, in the present manuscript we have first placed a proposition that in order to treat all higher order 
theories (including gravitation) on the same footing, one should fix acceleration instead of velocity at the end 
points. Once such proposition appears tenable, then neither Ostrogradski’s, nor Dirac’s technique can successfully 
handle the situation towards canonical formulation of Pais-Uhlenbeck fourth order oscillator action |4], which 
was presented by Mannheim and Davidson [23] . along with other higher order theories. This we demonstrate 
in section (2.1). On the contrary, incorporating such proposition in Horowitz’s technique it is found that the 
technique is unable to handle the counter terms in general. This has been demonstrated in section (2.2). In this 
context, we present a rigorous scheme for canonical formulation, by modifying Horowitz’s method in section (2.3). 
This technique treats gravitational and non-gravitational higher order theories on the same footing, removing 
the pathologies discussed. Although, the Hamiltonian found in the process is canonically related to the ones 
found using other techniques for the actions under consideration, in the appendix we have demonstrated that 
for more general action, the Hamiltonian may be quite different and is not related to the others under canonical 
transformation. Second, we demonstrate in section (3) that it is possible to regulate the pathology of branching, 
by associating higher order term in higher degree theory. 






2 Fourth order oscillator 


Since our aim is to explore the problem associated with higher order theory, meaning that the Lagrangian is a 
function of at least second derivative of a field variable in the form C — C{qi,qi,qi,qf,...) and the resulting Euler- 
Lagrange equation contains fourth or higher order terms, let us take up fourth order oscillator equation of motion, 
viz. 


q +{wl + wl)q + wiW 2 q = 0 , 

which is supposed to be found from the Pais-Uhlenbeck action [3] 


( 6 ) 


S'! = ^ J[f - {wl + wl)q^ + WiW2q^]dt. (7) 

Canonical formulation of the above action following Dirac’s constraint analysis was found by Mannheim and 
Davidson [23] as 

Hi = ^ + yPq + lj^{wl+wl)y‘^-'^WiW2q^, ( 8 ) 

where, y = q and Py is the canonically conjugate momentum. In fact, it is trivial to show that the same could 
have been obtained following Ostrogradski’s technique [5] also. However, under variation of the above action © 
one obtains 


<5S'i = 7 J [ 4 +{wl + wl)q + WiW2q]Sqdt + ^q6q, (9) 

using the condition (^qlendpoints = 0. Now, to perform canonical analysis, an additional degree of freedom is 
required, which is y = q. Therefore, to match boundary data with Cauchy initial data, yet another quantity is 
required to be fixed at the end points. Two options are now open. One is to set 5q\e.ndpoint = 0 and the other is 
to set 5q\e.ndpoint = 0. For the second choice, the action ([7|) should be supplemented with an additional term, so 
that it reads 


51 = ^ J[4^ - {wj + wl)q^ + wiW2q^]dt - ai 


= / Lidt — tJi 


( 10 ) 


where, tJi = ^qq. In Ostrogradski’s [^ (together with its generalization developed for non-singular gravitational 
Lagrangian by Buchbinder-Lyakovich and Karataeva |24] following Dirac’s constrained analysis) and of course 
in Dirac’s methods, first option is considered which apparently resolves the issue. But then, some problems 
appear. To better understand, let us consider the following oscillator action 


52 



^4^ + CU{q) 


dt, 


( 11 ) 


which, despite the presence of q term in the action, is not a higher order theory, since it does not produce fourth or 
higher order equation of motion. Clearly the first term is a total derivative term, which has been kept deliberately 
to administer the issue under discussion. Under variation (setting Sq = 0 at end points) one obtains 


SS 2 = J{Aq + CU')5qdt + ^ 9 ^(5qI endpoints, ( 12 ) 

which produces the equation of motion for ordinary harmonic oscillator (in the above, prime stands for derivative 
with respect to g), if one can regulate the end point data. With (5g = 0 at the end points, the endpoint data 




exceeds Cauchy initial data. So, the only option is to find a counter term for action (ED, so that the boundary 
term gets cancelled. The complete action then reads 


^2 = 


k Q .d. n ^ ^ 

^q^q-^q^+CU{q) 


dt — (72, 


(13) 


where, (T 2 = . Clearly, under integration by parts, q term gets cancelled with the counter term, and the action 

takes the form of usual oscillator action, despite the presence of second derivative term in the action. Thus, the 
need for a boundary term is realized, through the simple example cited above. Let us now modify action ED by 
introducing if term, so that it becomes the Pais-Uhlenbeck fourth order oscillator action with a total derivative 
term in the form 


^3 = 


n /c o ■ ■ .T o 

—(7^ + -q^ q - 

2 ^ 2^ 2 ^ 


CUiq) 


dt. 


Under variation one obtains 


(14) 


5S3 = 

J (B q + Aij + CU')Sqdt + ^ — f + Bq'^Sqlendpoints (15) 

{dd q ^())^(?|endpoints- 

Now if one fixes both q and q at the endpoints, boundary terms vanishes and there is no need for a counter 
term CT 2 . Appropriate field equations are retrieved and so apparently there is no problem. However, there is no 
reason as to why the counter term (T 2 disappears, when higher order theory is considered. Further, for H = 0, 
the counter term a 2 being absent, field equations are found only under the condition q and q fixed at the end 
points, and in the process, boundary data exceeds Cauchy data. lj. This is the same situation that one encounters 
in higher order theory of gravity, as mentioned in the introduction. Therefore, fixing both q and q at endpoints 
lacks mathematical rigor. 

On the contrary, if one chooses a counter term in the form = Bqq + and set Sq = 0 = Sij at the end 
points, then the complete action reads 


S. = 


B ..n k 2 ■■ A 2 y-iTTf \ 

+ q-_f + CU{q) 


dt — as- 


(16) 


Cauchy initial data now matches the boundary data and also, in the absence of higher order term (B = 0), the 
counter term for the harmonic oscillator action under present consideration ED, is retrieved. Hence we propose 
that one should fix acceleration and not the velocity at the end points for fourth order theory and in the process, 
gravity is treated on the same footing as other higher order theories. 


2.1 Canonical formnlation of Pais-Uhlenbeck action 

In view of preceding discussions, if the proposition of fixing acceleration at the endpoints instead of the velocity 
appears tenable, then obviously, Pais-Uhlenbeck action © should be replaced by (ED, which contains a counter 
term. The question is how shall we proceed for canonical formulation of the said action, handling the counter term 
appropriately? One option is to take derivative of ai and to enter it in the Lagrangian. In the process the action 
reads. 





2q q —{wl + wfjf + WiW 2 q^]dt. 


(17) 


^One would then resolve the issue of counter term appearing in lower order theory by adding a higher order term and then setting 
the higher order term to zero. In this way one can get rid of GHY term too. Clearly this is illegal. 









Under variation, action (O produces appropriate equation of motion setting Sq = 0 = Sq, at the end 
points. But, despite the fact that the above action produces fourth order Euler-Lagrange equation of motion, 
the appearance of third derivative term in the action functional requires yet another additional degree of freedom 
cc = q, in order to follow Ostrogradski’s technique [5] or to make Dirac constrained analysis [6]. In the process, 
Cauchy data exceeds boundary data. So both the techniques fail to handle fourth order oscillator action, under 
the current proposition. The other option is to take the action (US and to make use of an auxiliary variable 
Q, following Horowitz [7]. Although, the prescription was initially proposed for gravitational action, but since 
it is applied in minisuperspace model only, so it works for other higher order theories as well. Although, the 
prescription was initially proposed for gravitational action, but since it is applied in minisuperspace model only, so 
it works for other higher order theories as well. Although, velocity is kept fixed at the endpoints in this formalism 
too, it is indeed possible modify the technique by fixing acceleration instead, at the end points in Horowitz’s 
formalism. The prime issue of this scheme then reads: “first, one should choose an auxiliary variable varying 
the action with respect to the highest derivative of the (field) variables appearing in the action. Second, cast the 
action in canonical form, so that in the process of integration by parts, counter term gets cancelled. One can now 
vary the action both with respect to the fundamental variable q and an auxiliary variable Q. While variation of 
q yields appropriate Euler-Lagrange equation, variation with respect to Q returns its definition only. Einally, as 
the phase space structure of the Hamiltonian is found, the auxiliary variable Q should be replaced by the basic 
variable q, following appropriate canonical transformation”. The important point to note is that, q is treated as 
basic variable, but it doesn’t acquire the same status as q, since, neither it is required to vary the action with 
respect to g, nor it is fixed at the end points. With the above prescription, the auxiliary variable for action (HQD 
therefore is. 


Q 



in view of which the action may then be judiciously expressed as 


(18) 


Si 




(19) 


so that, under integration by parts, the boundary term oi gets cancelled and one is left with the following canonical 
action 






( 20 ) 


Q variation equation returns the definition of Q given in (1181) and q variation equation returns the fourth order 
oscillator equation of motion ([6|). The Hamiltonian reads 


Hi = -PqPQ + '^{wl + wDpI^ 


( 21 ) 


Now, phase-space structure of the Hamiltonian should be expressed in terms of the basic variables q and y = q 
instead of q and Q , as mentioned. This may be achieved by the following canonical transformation 


dLi dLi 


( 22 ) 


Hence the final form of the Hamiltonian is 


Hi = 


27 


VPq + - ^WlW 2 q‘^ 


(23) 


Note that this is the same Hamiltonian ((^51) obtained by Mannheim and Davidson following Dirac’s constrained 
analysis. In fact, all the three formalisms discussed here, viz. Ostrogradski’s [5], Dirac’s [6] and Horowitz’s [7], 
produce the same and correct phase space Hamiltonian for the Pais-Uhlenbeck oscillator action. However, while 
Ostrogradski’s and Dirac’s formalism keep velocity at the endpoints fixed, it is indeed possible to fix acceleration 
at the end points in Horowitz’s formalism. However, such a modified version of Horowitz’s formalism is found to 
be diseased from certain pathology, which we expatiate underneath. 






2.2 Pathology appearing in Horowitz’s formalism 

Although, Horowitz’s formalism produces correct phase-space Hamiltonian, keeping velocity fix at the end points, it 
is indeed possible to fix acceleration at the end points too. However, such modification often does not guarantee to 
account for the supplementary boundary terms correctly. Let us consider action ([T 6 )l to explore the situation. This 
action produces fourth order oscillator equation of motion provided, endpoint data are taken care of, appropriately. 
Plugging in the auxiliary variable following Horowitz’s prescription [7], 

Q = Bq+^e, (24) 

judiciously into the action (IT51) as. 


^3 = 


/ 





+ CUiq) 


dt 


- 0-3, 


(25) 


it may be cast in the following canonical form under integration by parts 


^3 = 



+ Qq- 0 - 3 . 



■ CU[q) 


dt 


(26) 


One can now trivially check that, Qq = Bqq+ ^q^q does not cancel the counter term <73. The situation worsen 
for a coupled action in the form 


S' = 7 


1-2 1 . .2 
12 '^! ■ 


CUiq,) + H<l2)q^qi 


+ ^92 - 1^(92) 


dt. 


(27) 


which is encountered in dilatonically coupled Gauss-Bonnet action associated with higher order term (say, ) in 
4-dimensional Robertson-Walker minisuperspace model. 


Finally, to apprehend the whole situation, let us cite a more concrete example in view of the following action 


S 4 = 7 


Iq^ + qnq -^Ae+CUiq) 


dt. 


(28) 


Under variation, one ends up with 


5S4 = -fJ [ 9 +Aq + CU']6qdt 
+ 7 [ 9"'^9 + 9^9 - QSq- AqSq] 


endpoints 


(29) 


Thus, the same fourth order oscillator equation of motion d®]) is reproduced, under the choice A = w, + W2, C = 
W 1 W 2 and C/(g) = , with appropriate end point data. However, following Ostrogradski’s or Dirac’s prescription 

(setting Sq — 0 — Sq at endpoints), one can’t retrieve the equation of motion, unless (5g|endpoints = 0. But then, 
endpoint data exceeds Cauchy data. Otherwise, one would require an additional degree of freedom to handle the 
situation. Nevertheless, in that case, Cauchy data exceeds endpoint dataH. On the contrary, the action may be 
supplemented by, 

^One might suggest that, set q and q at end points and supplement the action with a counter term , and thus the problem 

is resolved. But then the question arises, “at any stage if we require a counter term, then why shall we sacrifice the GBY term?” 













so that setting 6q = 0 = Sq at the end points, it reads, 


S'4 = 7 


/ 



( 31 ) 


/ 


L^dt — (J4. 


Now, following Horowitz’s prescription, the auxiliary variable Q is defined as 



(32) 


It now turns extremely difficult, if not impossible to cast the action (1311) in canonical form. As a result it is not 


known if the total derivative term cancels the supplementary boundary terms. 

2.3 A legitimate Canonical formulation 

At this end, we understand that, for higher order theory under consideration, if our proposition that Sq = 0 = Sq at 
the endpoints, appears tenable, then neither Ostrogradski’s nor Dirac’s canonical scheme works. On the contrary, 
although Horowitz’s formalism appears mathematically rigorous, it can’t handle the supplementary counter terms, 
in general. However, the pathology appearing in Horowitz’s formalism disappears if the action is integrated by 
parts, prior to the introduction of auxiliary variable. For example, the last action dMD under integration by parts 
yields 



(33) 


The above action (1331) takes the form of Pais-Ulhenbeck action (nni). Canonical formulation then follows as de¬ 
picted in section (2.1). The same technique resolves the problem of all higher order actions as well. 


In the case of gravity however, choice of appropriate variable is a precursor to handle the associated boundary 
terms [laiiniiiTiiiH]- As an example let us consider the following action 



(34) 


which, in the Robertson-Walker minisuperspace ([H reads 



(35) 


SttG 


3 



odd -b d^ -b ka) d^x. 


Under integration by parts one ends up with 



( 36 ) 



















Now introducing the auxiliary variable 


dC 

Q = — = 72f3aa, 
oa 


judiciously in the above action 


( 37 ) 


A,= 


+36/3 


3 

{—aa? + ka) + Qd — 


SttG 


144/3a 


(a^ + ky 


drxdt — 72/3 / add — I ax, 


it takes the following canonical form, under integration by parts 
714 = 


+36/3 


SttG 

(d^ + ky 


3 • 

{—ad^ + ka) — Qd — 


144/3a 


d^xdt + 48d^d^x. 


(38) 


(39) 


Nevertheless, one retains an additional boundary term. This is due to bad choice of co-ordinate, as we know, the 
basic variable is, hij = and not ‘a’. One can observe that under the choice hij = a^ = z, the action (1341) may 
be cast as 


A4 


3 




kz 

7^ 


+A:^i/z)] d^xdt 


IGttG 


yfzz 




which under integration by parts, becomes 


(40) 


714 = 


Z2 


SttG \ 4-y/i 

dt — 36/3 


+ fc-v/z ) + 36,0 


4Vi 


zz \ 


Now, plugging in the auxiliary variable 

Q = ^ = lSz^/z, 
oz 

the action reads 


714 = ^ 
+36/3 


SttG \ 4-y/z 

(K 


+ k-s/z I + Q 'z — 


Q^yfz 


36 




+ k^^z 


Z2 


dt — 36/3 


zz 


2^J ’ 


which under integration by parts yields the following canonical form 


7l4 = 


SttG \ Ay/z 


+ ky/z I — Qz — 


Q^yfz 

36 


+36/3 + k'^\Z^ 


dt. 


(41) 


(42) 


(43) 


(44) 


In the process, the problem with the boundary term disappears and one can now proceed to cast, what we call 
- the legitimate canonical formulation of higher-order theory, in which all the higher order theories (including 
gravitation) may be treated on the same footing. 



































3 Resolving the issue of branched Hamiltonian 

Since, we have been able to construct a viable canonical formulation, that fits both the modified fourth order 
oscillator action and gravity as well, let us turn our attention to resolve the issue of branched Hamiltonian. 
We have already mentioned that if a Lagrangian contains velocities with degree greater than two, the resulting 
Hamiltonian is multivalued, and the standard Legendre transformation to cast the action in canonical does not 
apply in general. Despite attempts, the problem associated with branched Hamiltonian has not been resolved 
uniquely. However, as already mentioned, such situation often appears with higher order theory. For example, in 
the context of gravity, higher order curvature invariants etc.) contain terms with higher degree in 

the field variables. We shall show that pathology of branching may be removed, if the action accompanies higher 
order term in addition, as presented in action ([^ . 


3.1 Generalized Pais-Uhlenbeck action 

Note that in the presence of an appropriate potential term, action ([5]) leads to Pais-Uhlenbeck oscillator action 
m accompanying a higher degree ( 9 '*) term. Therefore, for the sake of demonstration, let us take the modified 
version of the forth order oscillator action (EB, by associating with it a higher degree term, so that the Lagrangian 
C = C{qi^ (ji, (ji, qf ,..., qf) , and dub it as the generalized fourth order oscillator action. Such an action reads 


35=1 


lf+rq-^e + CU{q) + jq^ 


dt — a 


(45) 


In the absence of higher order (if' ) term, the above action suffers from the pathology of branching. However, in 
the presence of higher order term, integrating by parts yields 


85=1 

= / C^dt - 'yqij. 


lf-lAf+CUiq) + jf 


dt — 'Yqij 


(46) 


Introducing the auxiliary variable 


Q 



into the action as before and under integration by parts, one obtains 


(47) 


S. = 


-Qq-^- ^Acf + ^cu(q) + ^cf 


dt. 


The phase space structure of the Hamiltonian may then be obtained at-ease as, 

2 rmf \ . 

—pQ-jCU(q)-—pQ. 


H = -PQPq + §(; + '^Pq - l(^U(q) - 


(48) 


(49) 


Now replacing Q and pg by basic variables under canonical transformations, Q = Py and pq = —y, where, 
y = q, ss before, the Hamiltonian may be expressed in its final form as 

H = ^ + yPq + l[^y^-CU(q)-^y^). (50) 

Clearly, the velocities appear in the form of potential and the pathology appearing from branching is resolved. 


3.2 Higher order theory of gravity 

For the sake of completeness, let us turn our attention to higher order theory of gravity, yet again. Canonical 
formulation of gravity requires a 3-1-1 decomposition in 4-dimension, under which hij is treated as basic variable. 
In higher order theory of gravity, the extrinsic curvature tensor is treated as an additional variable, which 
essentially is the first time derivative of hij. Further, higher order theory of gravity requires to fix R at the 








boundary in addition, which in principle may be expressed in terms of hij and its second time derivative, viz. 
hij. To better understand, let us take Robertson-Walker minisuperspace model ([2]) , as an example. Apart from 
the 3-volume factor, since hij = , so instead of the scale factor a, one should treat z = as basic variable. 

Automatically, the additional variable is i = 2aa oc Kij. Therefore to reconcile boundary conditions, in addition 


to z, z should be kept hxed at the boundary. Now, 


R = 6( - + 4 + — I =3( - + 2- 


(51) 


Therefore, keeping z fixed is analogous to keep R fixed at the boundary, as already mentioned in the introduction 
0. Hence there is indeed no contradiction with higher order theory treated so far, since the choice of basic variables 
and boundary conditions reconcile. 


Next, it is required to see how multi-valued Hamiltonian appears in the theory of gravity, and may be handled 
taking into account higher order curvature invariant terms. It is well known that gauge invariant divergences make 
general theory of relativity non-renormalizable and a renormalizable theory of gravity in 4-dimension requires to 
modify Einstein-Hilbert action by incorporating curvature squared terms in the form [29j 

^5 = y 

However, analysis of linearized radiation reveals the presence of five massive spin-2 particles whose excitations 
are negative definite and therefore are ghosts. These ghosts destroy the unitarity and so no sensible physical 
interpretation of such a theory exists. In the absence of term ghosts disappear, but the theory becomes 

non-renormalizable. Although scalar curvature square term R^ is not responsible for the appearance of ghosts, 
yet it is a general practice to cast the gravitational action with a particular combination of higher order cur¬ 
vature invariant terms so that fourth derivative terms disappear from the field equations. Gauss-Bonnet term 
{R^ — AR/j^^R^’' + RapfiuR’^^^’^) was recalled for the purpose. As it is topological invariant in 4-dimensions, 
Lanczos-Lovelock gravity |16j was therefore constructed, which is realizable in dimensions higher than four. Nev¬ 
ertheless, it was found that the action contains higher degree term, leading to branched Hamiltonian and so the 
theory again lacks unitary time evolution of quantum states. This issue was resolved by modifying Lanczos- 
Lovelock gravity under the addition of a scalar curvature square term ( R?) |22] . In the following, we therefore 
take up yet another example in 4-dimension. 


R-2A 

2k 


+ PR^ + 


(52) 


The combination of higher order terms in the form TZ 2 = R — — R^) also produces second order 

field equations, although in Robertson-Walker mini-superspace ([2|) only. Hence, it might appear that the following 
action 


A6 = J yTgd^x + B^nl - - V{4>)^ 

~ J d^xy/hK, 


(53) 


should be free from all pathologies at least in Robertson-Walker minisuperspace metric and one does not 
have to bother about the boundary term d^xy/hKf'{R) along with additional boundary condition ((5i?| av = 0) 
appearing from f{R) theory of gravity, other than Gibbons-Hawking-York term [8]. Further, note that since in 
isotropic space-time, {Rf^^R^" — ^) is a total derivative term, so is usually replaced safely by R^ term. 

However, as mentioned, the theory is not renormalizable in the absence of R^^ term. But here, the presence of 
such a combination within the square root does not allow such replacement and hence both the terms are retained 
as such. Now, in terms of the basic variable (hij = a? = z), the Ricci scalar and the higher order curvature 

invariant term read R = 6(^ -I- y) and R 2 = 12 -|- ^ j respectively. After getting rid of the GHY boundary 


term under integration by parts, the action (1531) can therefore be expressed as 

-|- ky/z 

1 


A, = j 

_K V Ay/Z 


f z^ k \ 

144^2 

—5 - T ) 


- V{4>) 


3 


dt. 


(54) 


^In D dimension, R = {D — 1) + (D — 2)^j +6^ = 4 ^ ^ +6—^ , under the choice, z = a ^ . Particularly, in asynchronous 

frame ds^ = —dt^ + h^isdx^^dx^ , it is always possible to express R in terms of hij and hij . 








The presence of only first time derivative of the field variables clearly indicates that despite the presence of higher 
order curvature invariant terms, the field equations are only second order. However, the canonical momentum 


3 z 

Pz = — 7n= 


+ 144^2 




(55) 


involves z^ term, leading to the issue of branched Hamiltonian. Therefore, while perturbatively higher order terms 
destroy unitarity, nonperturbatively, higher degree terms destroy the same. However, such a worse situation may be 
resolved by adding yet another curvature invariant term where, 'R.\ = i?+-\/3(4i?^i,i?^'^ — W-) = 12 
to the action (15^ . so that the modified action reads 


M = 

J y/^d^x + B,nl + B2nl - 

~ ^{4’)^ + (T + s. 


(56) 


where, a is the GHY boundary term [5] and E = 4i3i J ^Ry/hKd^x, appears under variation, with an additional 
condition that 6R = 0 at the boundary. In fact, for B 2 = Bi, the above combination gives scalar curvature 
squared term R? = TZ\ + Ti,\ and the action simply reads 


A7 = 


'^R + b,r^ - - V{4>)] V^d^x 

tr + E. 


(57) 


A legitimate canonical formulation of the above action in the absence of the scalar field exists in the literature 
taking lapse function into account [mill]. Here we briefly discuss the outcome, in the context of the issue of 
boundary term, canonical quantization and the additional boundary condition in connection with Noether sym¬ 
metry. 


The issue of boundary term: 

In Robertson-Walker space-time dH) the action (l57l) , under the choice (B — 36Hi) reads 


At = 

3 / yfz'z 

K 


-|- ky/z ) -|- B 
dt -j- 


+ 


Ay/z yfz yfz 

(T -f El -|- E2. 


(58) 


In the above, the supplementary boundary term has been split as, E = 4i?i J(^R)Ky/hd^x = Ei -|- E 2 = 
4Hi j(^R)Ky/hd^x ABi f(^R — ^R)Ky/hd^x. Under integration by parts, the surface terms a and Ei get 
cancelled and the action takes the form, 


Afi = 


4yi 


-|- ky/z ) -|- H 


/ z^ kz^ 


= [ C^dt -|- E 2 , 


\4:y/z 2z^ 

dt E 2 


(59) 


which in view of the auxiliary variable 
^ dz 2 yi’ 


(60) 


1 n 














may be cast in the following canonical form 


Bkz^ 

+ 


2z^ 


4yi 

/I 15 

7in2^ 


+ ] — Qz — 




B 


-V 


(61) 


dt, 


and in the process, the boundary term E 2 gets cancelled with the total derivative term. Therefore, boundary 
terms have been taken care of appropriately, in two stages. The canonical momenta are 


3 i • 

PQ = -z, Pz = --Z - 1= ~ ^ 

2k y/z 


Bkz ■ 3 

= (jiZ'^ 


(62) 


Canonical quantization: 


The phase-space structure of the Hamiltonian is. 


AKy/z K 2 z2 2z2 


PzPQ 


Q^yfz Bk^ 


B 




+ V{(j})z^ = 0. 


(63) 


Now, in view of the definition of the auxiliary variable and its canonical conjugate momenta given in expressions 
dmi) and (I62p respectively, we are in a position to translate the Hamiltonian in terms of the basic variables 
{hijjKij). For the purpose, we replace Q by Py and pq by — y, where y = z. The Hamiltonian therefore takes 
the form 


H = yp^ + ^pI 
Bky"^ 


+ Vz^ — 


^ pI 3 y^ 3 

1- 3 —I— -— —ky/ Z 

2zi 4 k y/z K 

Bk^ 


2zi 




(64) 


which is constrained to vanish. Canonical quantization then leads to the following Schrodinger-like equation 
ih 94' dA/ 



4- - 


924 - 


VW - ^ 

y 2r] 


2r]ty dcjP' 

Bky Bk'^ 




4', 


(65) 


where, y = zi = - the proper volume, acts as internal time parameter and n is the operator ordering index. 

Clearly one ends up with a hermitian effective Hamiltonian operator with unitary time evolution of quantum 
states. In the process, the issue of branched Hamiltonian has also been resolved. It has been shown [23 [28] that 
the above quantum description leads to well-posed probabilistic interpretation, fixing the operator ordering index 
to n = — 1. Further, the semiclassical wave-function has been found to be oscillatory, indicating that the region 
is classically allowed and is strongly peaked about a set of inflationary solutions admissible by the classical field 
equations. 


Noether symmetry: 

First of all let us mention that in the process of canonical formulation of higher order theory following Os- 
trogradsi’s and Dirac’s techniques, it is not possible to explore Noether symmetry at least in the Lagrangian 
formalism. However, this may be expatiated under the present technique of canonical formulation. Noether sym¬ 
metry for finding the form of the scalar potential corresponding to higher order theory of gravity in addition to 
non-minimally coupled scalar-tensor action already exists in the literature [80] . It was there observed that the clas¬ 
sical solutions of the field equations are obtained directly from the Noether equations. For the sake of simplicity. 



















here we ignore the Einstein-Hilbert part of the action (EZl, which is of-course ignorable in the very early universe. 
In the Noether symmetry approach, the lift vector X acts as an infinitesimal generator of Noether symmetry in 
the tangent space (z, i, Q, Q, (j), (j )), and is introduced as follows, 


X = a 


dz 


dQ 


d . d . d ■ d 

■ 7 — + aTTT + +a -^ + P - r. 

dz dQ dcj) 


d(j) 


( 66 ) 


The existence condition for symmetry, £xL = XL = 0, then leads to the following master Noether equation, in 
view of the canonical action dSI]), 

2,BkP BP 3/02 X ‘ 

- “^Qy/z-pz^V' + {a,zZ + a,QQ + a^)(^^^ ( 57 ) 

- Q) - {l3,zZ + /3,qQ + I3'4)z + 

+ y'0)0z5 = 0. 


Equating coefficients {z^,Q'^,(j£,zQ,Q(j),(j)z) and the terms independent of time derivative to zero, as usual, we 
obtain following set of partial differential equations, viz., 


3 5 

Bk(^za^z — = 0 = a^g = 3a + 4zy' 

ct,z + /3.Q = 0 = z^y,g - a' = Bka' + zip' - z^p^z 


( 68 ) 


Q2 sfc2 3 ; 


+ %Q^z£zipV' = ^. 

JD 


These set of equations may be solved for fc = 0 to yield. 


41 41 

a = -yz; P = —Q; p = Ip 


V{P)=VoP^ + —; Q = 


^^/B 


(69) 


^z, 


where, I and m are the separation constants. In view of the definition of the auxiliary variable (l60)) . R{= ^=) 
turns out to be a constant, leading to de-Sitter or anti-de Sitter solution depending on the separation constant m. 
The conserved current 


41 ■ 41 -3 

X = —Qz - —Qz + Ippz^, 
o o 


(70) 


satisfies the field equations and everything is well behaved. Since, Noether symmetry demands R to be constant, 
it is trivially constant at the boundary. Thus, in view of Noether symmetry, R^ gravity admits de-Sitter or anti 
de-Sitter solutions and nothing else. In a recent article |3I] . it was shown that Noether symmetry of higher order 
theory of gravity administers the same condition {R = constant) both in isotropic and anisotropic space-times. 
Thus, in view of Noether symmetry, we get an answer to the question “why higher order theory of gravity demands 
R should remain fixed at the boundary?”. 


4 Concluding remarks 

As mentioned in the introduction, there are several reasons which favour to keep the Ricci scalar (R) fixed at 
the boundary, for higher order theory of gravity, in addition to the metric {g^v)- Of particular importance is, 
scalar-tensor equivalence of higher order theory of gravity requires to fix R at the boundary. This means keeping 
the acceleration fixed at the end points, in disguise. On the contrary, Ostrogradski’s or Dirac’s techniques for 
canonical formulation of fourth order theory, tacitly assumes velocity should be fixed, instead. In the process, not 
only GHY term, which is related to the entropy of a black hole gets lost, but also there exists a clear contradiction 













in handling different higher order theories. 


To remove such contradiction and to treat different higher order theories on the same footing, we propose that 
acceleration, in addition to the coordinate should remain fixed at the endpoints. Such proposition, as mentioned, 
however is not compatible with Ostrogradski’s and Dirac’s techniques. Nevertheless, Horowitz’s technique of 
associating an auxiliary variable, although originally developed for canonical formulation of higher order theory 
of gravity, appears to work in general with the said proposition. Unfortunately, with such proposition, Horowitz’s 
formalism also often fails to handle supplementary boundary terms. In this connection, we have modified the 
technique and in the process, a one-to-one correspondence between canonical formulation of higher order oscillator 
action and gravity has been established. Further, in the appendix we have shown that the phase-space Hamiltonian 
constructed following such technique is different and is not related to the others under canonical transformation. 
The modified version is nnique in this sense. 

The presence of higher degree term in velocity leads to branched Hamiltonian and tells upon unitary time 
evolution of quantum states. Despite attempts, there does not exist a unique technique to resolve the issue. 
Under such circumstances, we have also demonstrated the fact that, the pathology due to the appearance of 
such multivalued Hamiltonian may be regulated if it appears with higher order theory. The technique was earlier 
applied to alleviate the pathology of branched Hamiltonian appearing in Lanczos-Lovelock gravity, by adding a 
scalar curvature squared term [22) . Here, the technique has been applied in a theory of gravitation, which appears 
to be free from all pathologies in connection with renomalizability and the appearance of ghost degree of freedom, 
of-course in isotropic and homogeneous cosmology. However, the problem has been resolved by adding a higher 
order term. The technique is well posed, leading to Schrodinger like equation with straight forward probabilistic 
interpretation, which is lacking in Wheeler-DeWitt equation. Semiclassical approximation had been found to be 
strongly peaked about classical inflationary solution mm- 

In the process of canonical formulation following Ostrogradski’s or Dirac’s technique, the Lagrangian is left 
alone, and indeed there is no way to explore Noether symmetry at least in the Lagrangian formulation. Use of 
auxiliary variable on the contrary, is well suited for the purpose. Now, one way to resolve the issue of dark energy 
is to modify the left hand side of Einstein’s equation by adding higher order curvature invariant terms. Usually, 
one chooses F{R) by hand or following the reconstruction scheme. Recently, in an attempt to find a form of 
F(R) by invoking Noether symmetry, both in isotropic and anisotropic space-time, it was observed that Noether 
symmetry demands R must remain constant, resulting in de-Sitter/anti de-Sitter solutions [31]. This at least 
gives a signal that nothing other than R could be fixed at the boundary. Here again we observe that Noether 
symmetry requires R — Rq - a constant, and as such administers the additional boundary condition dRjav = 0. 
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A Different phase-space Hamiltonian 

In the appendix, we cite an example to demonstrate that not all phase-space Hamiltonian corresponding to higher 
order theory, constructed out of different canonical formalisms are related through canonical transformation. Let 
us take the following higher order action. 



(71) 


Now, choosing d = x, the above Lagrangian 


L — {x + aY f{4>) + + 


,2 


2 


2 




(72) 


1 n 





becomes singular and so instead of Ostdtrogradski’s technique, one should follow Dirac’s constrained analysis, 
fixing a and d = a; at the boundary. Introducing Lagrange multiplier A, one then writes 


L — f(x + a)^ ~2 ~2 -^ ~ 

Canonical momenta are. 

Pa; = 2f{x + o). Pa = A, p 0 = 0, Pa = 0. (74) 

The constrained Hamiltonian then reads 

2^2 
He = ^ — apx - ^ 

Clearly we require two primary constraints involving Lagrange multiplier or its conjugate viz. 


0i = Pa - A, 02 = PA- (76) 

The primary Hamiltonian therefore reads 

2 ^2 

Hpi = ^ - apo;+ ^ - ^ + Xx + Ui{pa- X) + U 2 P\. (77) 

In the above ui and U 2 are Lagrange multipliers and the Poisson bracket {a,pa} = {a;,pa;} = {A,pa} = 1, hold. 
Now constraint should remain preserved in time. Therefore 


2 

01 = {01, Hpi} = -U2 +Pa; + ^ 0i{0i, Ui} = 0 

2 

02 = {02, Hpi} = (ui - x)px + ^ 0i{02, Ui} = 0 

i=l 

Since, Constraints must vanish weakly in the sense of Dirac, therefore ui 
Hamiltonian therefore reads 


(78) 


X and U 2 = Px ■ The modified primary 


Hp2 = xpa - apx + ^ + ip^ - + V + PxP\. 


(79) 


Now since 02 = { 02 , Hp 2 } = 0, so pA = 0. Hence, the phase-space Hamiltonian finally takes the form, 


TJ I Px I 

Hd = XPa - aPa, + ^ + 


I 


pI 


V. 


(80) 


The same phase-space Hamiltonian is produced following Horowitz formalism also. However, the so called legiti¬ 
mate canonical formulation lead to 


HL=xpa + ^ + ^p| + 2f'axp^ - ^ [(1 - 2f'af - 4f]x^ - fa^+ V (81) 

While, Hamiltonian (1501) does not require, the fourth term in the above expression (ISTl) . requires operator ordering, 
which is different for different forms of /(</>). Clearly, therefore. Ho and H^ are not related through canonical 
transformation. Such situation is encountered due to non-minimal coupling which appears in the case of any 
non-minimally coupled higher order theory of gravity. 


